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Abstract 

Yang-Baxter system related to quantum doubles is introduced and 
large class of both continuous and discrete symmetries of the solution 
manifold are determined. Strategy for solution of the system based on 
the symmetries is suggested and accomplished in the dimension two. The 
complete list of invertible solutions of the system is presented. 



1 Introduction 

The Yang-Baxter equations, both constant and spectral dependent proved to 
be an important tool for various branches of theoretical physics. They represent 
a system of cubic equations for elements of x N'^ matrix R and can be 
written in the well known form 

R12R13R23 — -R23-Rl3-Rl2- (1) 

Even though many solutions are known for all types of the Yang-Baxter equa- 
tions ^, 1^, ^, ||] , until now the complete solution is known only for the constant 
Yang-Baxter equations in the dimension two i.e. matrices 4x4. 

Various extensions of the Yang-Baxter equations for several matrices, called 
Yang-Baxter systems, appeared in literature 0, ^. The constant systems 
are used mainly for construction of special Hopf algebras while the spectral 
dependent solutions are applied in quantum integrable models. Examples of 
both types together with their particular solutions were presented in pO| . 

As the Yang-Baxter systems usually contain several Yang-Baxter-type equa- 
tions, it is convenient to introduce the following notation; The Yang-Baxter 
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commutator [R, S, T] of (constant) x matrices R, S, T is N"^ x matrix 

[R, S, T] :— R12S13T23 — T23Si3Ri2. (2) 
In this notation, the constant Yang-Baxter equation is written as 

[R, R, R] = 0, (3) 
while e.g. the system for quantized braided groups |^ reads 

[Q,Q,Q] = 0, = 0, (4) 

[Q, R, R] = 0, [i?, R, Q] = 0. (5) 

The complete set of invertible solutions of this system in the dimension two is 
given in . 

The goal of the present paper is to give a complete solution of another, more 
complicated constant system related to quantum doubles i.e. special quasitri- 
angular Hopf algebras constructed from the tensor product of Hopf algebras by 
defining a pairing between them. 

2 WXZ system and its symmetries 

In the paper a method of obtaining the quantum doubles for pairs of FRT 
quantum groups is presented: Let two quantum groups are given by relations 

m 

W12U1U2 = U2U1W12 
Z12T1T2 — T2T1Z12 

where W and Z are matrices N'^ x iV^ satisfying the Yang-Baxter equations 

[W, W, W] = 0, (6) 
[Z,Z,Z]^0, (7) 

and suppose that there is a matrix X that satisfies the equations 

[W,X,X]^0, (8) 
[X,X,Z]=0, (9) 

Then the relations 

X12U1T2 = T2U1X12 (10) 
define quantum double with the pairing 

<f/l,T2 >=Xi2. (11) 

On the other hand, the equations (H)-(|^), that in the following we shall call 
the WXZ system, are constant version of the spectral dependent Yang-Baxter 
systems for nonultralocal models presented in |p^ . 

There are two natural questions related to the WXZ system: 



2 



• Is there a matrix X such that for any pair of matrices W, Z that solve the 
Yang-Baxter equations the triple {W, X, Z) solves the WXZ system? 

• Is there for any matrix X a pair of matrices W, Z such that the triple 
{W, X, Z) solves the WXZ system? 

Answers to both these questions are positive because the following two simple 
propositions hold 

• Let W, Z are arbitrary solutions of the Yang-Baxter equations [W, W, W] = 
0, [Z, Z, Z] = 0. Then the triple {W, X = 1,Z) is a solution of the system 

• Let X is an arbitrary matrix N'^ x iV^ and P is the permutation matrix. 
Then the triple {W ^ P,X,Z = P) is a solution of the system 

Other simple solution of the WXZ system is {W ^ R, X = R, Z ^ R) , where R 
is an arbitrary solution of the Yang-Baxter equation. 

Besides the above mentioned, one can find solutions of the WXZ system from 
the knowledge of solutions of the system Namely, if (Q, R) is a solution 

of the system of Yang-Baxter type equations then {W = Q, X — R, Z = 

Q) and {W = Q,X ^ R,Z ^ R+QR-) are solutions of the system (|)-(|). 

Solution of the system is essentially facilitated by knowledge of its sym- 
metries. It is easy to check that the set of solutions is invariant under both 
continuous transformations 

W = Lu{T®T)W{T®T)-^ 

X' = £,{T(g)S)X{T(g)Sy^ (12) 
Z' = C{S^S)Z{S^S)-^ 

where 

o.,^,CeC, T,Segc{N,C), 

and discrete transformations 

{W',X',Z') = {W^,X^,Z^) (13) 

iW',X',Z') = {W'\X,Z''), a = id,#, b = id,# (14) 

{W\X',Z') = {W^,X-,Z^), c = +,-,d=+,- (15) 

{W, X', Z') = {Z", X+, W^), c = +, -, d = +,- (16) 

where is transpose of Y , y+ := PYP, P being the permutation matrix, 
Y- := Y-\ Y^"^ := Y, Y* := {Y+)- = {Y-) + . 
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These symmetries can be of course composed, so that one can define e.g. the 
" antidiagonal transposition" 

F"* := (cr® cr)r^(CT(g)a) where (T = J (17) 

and the symmetry 

{W, X', Z') = (W*, Z"*) (18) 
Beside that there are useful " conditional symmetries" that can be expressed 

as 

Lemma 1: Let W,X,Z are matrices x that solve WXZ system 
(||)-(||) and T, A are matrices N x N such that 

[W, T(g)T]=0, (19) 

XiT ® 1) = (A 1)X or (T 1)X = X{A 1). (20) 



Then 



and 



W ^W, X' = {T® 1)X, Z' = Z (21) 
W ^W, X' ^ X{T (g) 1), Z' ^ Z (22) 



W' ^{T® l)WiT X' = X, Z' = Z (23) 

solve the WXZ system as well. 

Lemma 2: Let W,X,Z are matrices x that solve WXZ system 
(|)-(|) and S, A are matrices N x N such that 

[Z, S(g>S]^0, (24) 

X{S (g) 1) = (A (g) 1)X or (S" (g) 1)X = X(A g) 1). (25) 



Then 



and 



W' = X' = {S(E) 1)X, Z' = Z (26) 
W ^W, X' = X{S (g 1), Z' = Z (27) 

W' ^W, X' ^ X, Z' ^ {S®l)Z{S®iy'^ (28) 



solve the WXZ system as well. 

Proofs of these lemmas can be done by direct calculations. There are other 
symmetries of the WXZ system that are extensions of the twisting transforma- 
tions of the solutions Yang-Baxter equation but we are not going to use them in 
the following as it seems that they do not produce equivalent quantum doubles. 
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3 Solution of the WXZ system in the dimension 
two 



Even for the lowest nontrivial dimension two, solution of the system (^)-(^ rep- 
resents a tremendous task, namely solving 256 cubic equations for 48 unknowns. 
That's why it is understandable that the assistance of computer programs for 
symbolic calculations is essential in the following. On the other hand it does 
not mean that one can find the solutions by pure brute force, namely applying 
a procedure SOLVE to the system of the 256 equations. Moreover, the interpre- 
tation of the result would be extremely difficult as hundreds of solutions would 
appear, many of them related by the symmetries. 

There are several strategies for solution of the problem given above. All 
of them are based on the symmetries of the system and the knowledge of the 
complete set of solutions of the Yang-Baxter equation in the dimension two ^ . 

One possible (and obvious) strategy is solving the equations for all 

pairs of matrices W, Z in the Hietarinta's list of solutions of the Yang-Baxter 
equation. By this way we reduce the problem to 128 quadratic equations for 
16-22 unknowns (depending on the number of parameters in the solutions of the 
Yang-Baxter equation). 

Another strategy is to use the symmetry ( p^ to simplify the matrix X as 
much as possible, then solve the (linear in W or Z) equations (|8|),(^) and finally 
solve the Yang-Baxter equations (|^),(0) for W and Z. This is an analogue of 
the strategy accepted in [|ll| for solving the system (||)-(||). 

Both the mentioned strategies yield sets of equations and unknowns that are 
still too large to be solved by the computer program. The strategy that seems 
to work is a combination of the two previous. 

First we find the matrices X that solve the equation (||) for each invertible 
W that belongs to the Hietarinta's list. In the next step we solve the equation 

® ■ n . 

and finally determine Z from ([7|) using the results of the previous step. As 
the list of solutions is quite large it is essential in each step to factorize its results 
by the symmetries (|l2|)-(|l6|). 

Detailed description of this procedure and results are presented in the fol- 
lowing subsections. 

3.1 Solutions of the Yang— Baxter equations 

As it was mentioned, the important point is the knowledge of all invertible 
solutions of the constant Yang-Baxter equation [W, W, W] = 0. Their list up 
to the symmetries, denoted S, consists of twelve (intersecting) subsets |2[ || 
parametrized by up to three complex numbers. 

S := {R3.i{r, s, i), i?2.i(r, s), i?2.2(?', s), i?2.3(u, v),Ri,iir), i?i.2(s), Ri.si^), 

^Reduce 3.6 and Maple V was used 



5 



RiA{t), Ro.i,Ro.2, Ro.3, RoA \r,s,teC\ {0}, u,v eC} 



(29) 



The matrices R are defined in the Appendix 5.1. 

For later use we shall define three subsets of S, namely 



S5V := {i?2.i(r,s),i?2.2(r,s),i?3.i(r,s,t)| r, s, i G C \ {0}}, (30) 

Ssv ■= {R3.i{r, s, i), i?2.i(r, s),R2.2{r, s), -Ri.i(r), i?i.2(s), RiA{t), 

i?o.i, i?o.2, i?o.3, Roa\ r, s, t e C \ {0} } (31) 
SsT := {i?3.i(l, 1, 1), i?2.3(M, v), Ri.siu), i?o.i| u,ve C}. (32) 
They are the sets of solutions of the " five-vertex" form 



R: 



'5V 



"eight -vertex" form 
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■sv 



and "special triangular" form 
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(33) 



(34) 



(35) 



It is easy to check that 



S — Ssv U SsTi Szv C Sg,v 



(36) 



The reason for this partition is that these subsets have common symmetries and 
common matrices X that solve the equation (||). Namely, each element of S^v 
is invariant under 

R' ^ {A® A)R[A® A)-^ (37) 

where 

a 
h 



A 



(38) 
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while each element of Ssv is invariant under (^7|) where 



n 1 ) (39) 



-a 

and each element of Sst is invariant under (^7|) where 

--ill)' w 

These invariances imply the existence of solutions X = A(E) B oi the equation 
(H) common to all elements of the subsets S^v, Ssv or Sst- Actually we shall 
find that slightly more general solutions are generic for each of the subsets. 

3.2 Solutions of the WXX equation 

The goal of this subsection is to find all invertible matrices X that solve the 
equation 

[W, X,X]=0 (41) 

for all invertible solutions W of the Yang-Baxter equation (^ . First of all we 
can exploit the symmetry (O) to reduce the matrix X to one of the forms 



Ai, . . . , Ai4 given in the Appendix 5.2, Then we choose a solution W of the 
Yang-Baxter equation and solve X from the equation ( ^) . The list of the 



solutions up to the WXZ symmetries (M), (M),(N), (22) can be described 



in terms of 23 classes of matrices parametrized by complex numbers and given 



in the Appendix 5.3 



As mentioned in the previous subsection, there are solutions common to all 
elements in subsets S^v, Ssv or Sst- Namely: 

• liW e Sqv then X^{a, b, c) and ^(a, b, c) solve (^) for aU a,b,ce C. 



• If e Ssv then A^(a, ei, £20), and A||(ei, a, eio) solve (41 ) for all a e C 
and ei, £2 = ±1- 

• If e Sst then A|(a, 6, c) and X^{a, b, c) solve (|l|) for aU a,b,ce C. 

We shall call these solutions generic. Notation of the X-matrices above as well 
as below corresponds to that in the Appendix |5.3| . 

Below we display all invertible non-generic solutions of (|4l]). The list is 
complete up to the WXZ symmetries (Hj), (|l|),(^, (|2|), (||). That's why we 
also present the transformations that leave given W invariant or form invariant, 
i.e. invariant up to the change of parameters. 

The ranges of parameters are 

r,s, teC\{0}, u,v,a,b,c,d,e, f,g,h e C, e,ei,e2 = ±l. 



List of non generic solutions of the equation (41) 
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. W^R3.i{r,s,t). 

This matrix is invariant w.r.t. (|3^),(^8[) and 

i?3.i(r,s,tf = i?3.i(r-,s,t), R3.i{r,s,t)-^ = R3.i{r-\s-\t-^), (42) 

i?3.i(r,s,t)+ =i?3.i(s,r,t) (43) 

(T®r)i?3.i(r,s,t)(T®T)-i=ti?3.i(^,^,^), where T= ( J J )■ 

(44) 

The non-generic solutions of the equation [W, X,X] =0 up to the sym- 
metries exist only for special values of the parameters r, s, t. 

1. VK = i?3.i(s,s,l), X = X^{a,b,c) 

2. VK = i?3.i(s,s,l), X = jqia,b,c) 

3. VK = i?3.i(s,s,l), X = ^(a,6) 

4. VK = i?3.i(s,-s,l), ^ = ^(a) 

5. = i?3.i(l/s,s,e), X = Jf|(s, a, esa) 

6. VK = i?3.i(e,e,-e), X ^ X^{e,b, -b, c) 

7. VK = i?3.i(-l,-l,l), X = JCg(a, 6, c, d) 

8. VK = i?3.i(l,l,l), X = X^ia,b,c,d,e,f,g,h) 

9. W^ = i?3.i(l,l,l), X = ^(a,5,c,d,e,/,5,/i) 

• W^R2.i{r,s) 

This matrix is invariant w.r.t. (|3^),(^8[) and 

R^\{r, s) = i?2.i(s, r), i?2.i(r, s)-^ = R2.i{r'\ s'^) (45) 

{T <E>T)R+,{r,s){T (g,Ty^ = R2.i{r,s), where T = ( J J )■ (46) 

The non-generic solutions of the equation [W, X,X] =0 up to the WXZ 
symmetries are 

1. W = i?2.i(r, s), X - X^ir-\a, sa) 

2. = i?2.i(s, s), X = JCg(s-i, a, sa, b) 

3. VK = i?2.i(l,-l), X-4i(a,6) 

4. W^ = i?2.i(i,i), X- J^g(a) 

5. W — i?2.i(i,— i), -'^ = -''^(a) (coincide with the case 4 in 14^ = 

i?3.l)- 

6. W = i?2.i(e, e), coincide with the cases 1-3,5,7-9 in = i?3.i. 
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• W = R2.2ir,s). 

This matrix is invariant w.r.t. (|3^),(^8[) and 

^2*2(^,5) = -r.si?2.2(-r-\-s-i), R2.2ir,s)-' ^ R2.2{r'\s-^) (47) 

{T^T)R+^{r,s){T<»T)-^ ^ ~rsR2.2hs-\-r-^), where T = J J ^ . 

(48) 

The non-generic solutions of the equation [W, X,X] =0 up to the WXZ 
symmetries are 

1. W = i?2.2(r, s), X = X^ir-\a, ~ar-^) 

2. VK = i?2.2(s,s), X = X^{s~\a,-s-^a,g) 

3. = i?2.2(l, —1), R2.2{i, i) coincide with cases 3,4 in = i?2.i- 

• W^R2Au,v) 

This matrix is invariant w.r.t. (^), ( ^ ) and 



i?5*3(u,«) =i?2.3(u,«) (49) 

(T T)R2,l{u, v){T ® T)"^ = i?2.3(u, 2m - w), where T = 



(T T)R+^{u, v){T (g)T)-^ ^ i?2.3(l/w, w/w^), where T : 



1 

-1 

(50) 

u 
1 

(51) 

and there are no non-generic solutions of the equation [W, X,X] =0 up 
to the WXZ symmetries. 

• W = Ri.iir). 

This matrix is invariant w.r.t. (|3^),(^9|) and 

Rl^{r) = Ri,i{r), i?i.i(r)+ = i?i.i(r) (52) 

(T^r)i?^J(r)(T®r)-i = -ii?i.i(-r), where T= ( J • )■ (53) 

The non-generic solutions of the equation [W, X,X] — Q up to the WXZ 
symmetries are 

1. i?i.i(r), X = Xg(r, -e, er, e) 

2. - i?i.i(r), X = X^ia, e, i±i), r ^ 0, ±1 

3. VF = i?i.i(i), X = jqi](a) 
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This matrix is invariant w.r.t. (p^), (p9|) and 



(T ® T)R~l{s){T (g, r)-i = i?i.2(l/s), where T = J ^ ^ 



(54) 



and there are the following non-generic solutions of the equation [W, X, X] — 
up to the symmetries: 



-1 ^) 



1. W^Ri.2is), X = X|(l,e,-e) 

2. M/ = i?i.2(s), X-^(s-i,e,-es- 

3. W = i?i.2(s), X = ^(a, b), a^-P = s + l 
W = Ri.3{u) 

This matrix is invariant w.r.t. (|3^), ( pO| ) and 



i?+3(u) = i?j;;J(u) = (T T)i?i.3(u)(r r)"\ where T 



1 
-1 

(55) 

and there are the following non-generic solutions of the equation [W, X, X] - 
up to the symmetries: 

1. W^Ri.siu), X = X^{a,b,b-u-l,-ua), u ^ -1 

2. W^ = i?i.3(l), X^X^ia,b) 

W^RiAit). 

This matrix is invariant w.r.t. (|37|),(p8|) and 



RUit) = Ri.iit), Rli = Ri.iit), Ri.lit) = RiAit-'), (56) 

(T ® T)RiA{t)(T (g> T)-^ = RiA{t), where 7" = J J ^ • (57) 

The non-generic solutions of the equation [W, X,X] — up to the sym- 
metries are 

1. W^RiAit), X^X^{a,e,a) 

2. W ^RiAit), X = X^ia,l,-a) 

3. RiAit), X = X^{a) 

4. RiAit), X = X^{a) 

5. = i?i.4(±l). This matrices can be transformed to i?3.i(—l, —1, 1) 
by symmetry transformations so that these solutions are equivalent 
to the cases 1,2,3,5,7,8 (s -1, e -1) of W ^ Rs.i- 
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W = Ro.i 

This matrix is invariant w.r.t. (|3^), ( ^ ) and 

^0*1 = -Ro.i, -Rai — Ro.i, (58) 

(r(^r)i?(^}(r®T)-i = i?o.i, whereT= J (59) 

and there are the following non-generic solutions of the equation [W, X, X] 
up to the symmetries: 

1. W^Ro.i, X ^X^{a,b,c,e) 

2. W = Ro.i, X = Xa{a,b,c) 



W = Ro.2 

This matrix is invariant w.r.t. (|3^), ( ^ ) and 

R'a.2 — Ro.2^ Rq,2 = Ro.2^ (60) 

(T T)Rq I{T r)-i = i?o.2, where T = J ^ ^ (61) 

and there are the following non-generic solutions of the equation [W, X, X] 
up to the symmetries: 

1. W^Ro.2, X = X%{-l,e,-e) 

2. W = Ro.2, X = J(0(e, 6, 0, 1) 

W = Ro.3 

This matrix is invariant w.r.t. (^), ( ^ ) and 
Rl^ = Ro,3, (T®T)i?+3(r®T)-i =i?o.3, where ( 1 )' ^"^^^ 

2(T (g) T)Rq I{T (g) T)-i = i?o.3, where ^ = J ^ ^ (63) 

and there are the following non-generic solutions of the equation [W, X, X] 
up to the symmetries: 

1. W^Ro.3, X = X^{a,-l,a,i) 

2. W^Ro.3, X = Xaia,b) 



W = RoA = P 

Arbitrary matrix X solves the equation X,X]~Q. 
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3.3 Solution of XXZ and ZZZ equations for given X 



In the preceding subsection we have found all pairs {W,X) of invertible 4x4 
matrices that solve the equations (^), (||). The solutions are given in terms of 
the matrices X presented in the Appendix 5.3. The last step to do for solving 
the WXZ system is to find the invertible 4x4 matrices Z that for given X solve 
the equations 

[X, X, Z] = 0, [Z, Z, Z] = 0. (64) 

Trivial solution for any X \s Z = Rqa = P = permutation matrix. Below we 
display the list of nontrivial solutions of the equations (p4[) for all X-matrices 



from the subsection 5.3 



Let us note that in spite of the fact that the matrices Z must solve the Yang- 
Baxter equation, they need not belong to the list of solutions in the Appendix 

because in general we have not at disposal all the transformations ( |69[ ) up 
to which that list is complete, because they might already be used as a part 
of WXZ symmetries @, (|l4|), (|l|) to fix the matrices W and X. The only 
symmetry that we can use in general is 



(W',X',Z') = {W,X, Z*). 



(65) 



Beside this we can use the WXZ symmetries that leave X (and corresponding 
W) invariant. For some matrices X the whole sets of the solutions of the Yang- 
Baxter equation S or S^v or S^v or Sst solve the equations (|6j). Due to the 
restricted symmetries of matrices Z 

it is convenient for classification of solutions to introduce another special set 
of Yang-Baxter solutions 

TIbv := {-R3.i('', s, t),R2.i[r, s),R2.2{r, s), Qi?i.i(r)Q~\ 

QRi.2{s)Q-\QRi.2{s)^Q-\QRiA{t)Q-\ 
QRo.iQ ^,QRo.2Q ^iQRo.sQ ^,QRq.2Q ^tQRq.sQ 

/ 1 p/q^ \ 
1 p 
p 1 



V pq 







/ 1 p/q^ \ 
-1 p 
p -1 



1 / 



\ pq 







(66) 



1 / 



where 



q,r,s,t eC\{0}, pe C\{1,-1} } 



Q = D{q) = diag{l,q,q,q^) 



(67) 

This is the set of the " eight- vertex" solutions up to the symmetries R' — R'^ . 
The last two matrices can be transformed to R^.i{j^tj^A) a-nd Ria{jz^) 
by similarity transformations. 
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Below we display nontrivial invertible solutions of (JG^) . The list is complete 
up to the symmetries (p^-(p^, (|2^), (p7|). The ranges of parameters are 

q,r,s,t,x e C\{0}, u,v,a,b,c, f, g, h e C, e,ei,e2=±l, 

and S, Ssv, S^v, Sst are defined by (p9|)-(|3^). 

List of nontrivial solutions of the equation (|64|): 

• X = jq{a,b,c) 

The nontrivial solutions for general values of a,b,c are 

1. ZeS^v 

For special values of a, 6, c we get the following solutions 

2. (a = -1 A c = ±6) or (a = 1 A c = -6) : Z e TZsv 

3. a=lAc = fe: Z eSU {Ri,2{s)}- 

The nontrivial solutions for general values of a, 6, c are 
1. Z ^ Sst 

There are no other solutions for special values of a, 6, c. 

. X = J(|(a,6,c) 

The nontrivial solutions for general values of a,b,c are 

1. ZeS5v 

For special values of a, 6 we get the following solutions 

2. a -1 A -c : Z £ TZsv 

3. a=lAc = 6: Z eSU {R[,2{s)} 

• X = XJ^{a,b,c) 

The nontrivial solutions for general values of a,b,c are 

1. Z (z Sst 

There are no other solutions for special values of a, b, c. 

• X = X^{a,b,c) 

There are no nontrivial solutions for general values of a, 6, c. For special 
values we get 

1. 6 = -l,a = 1/c : Z = Ri.i{i) 
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2. b= -l,a=l/c : Z = Ria{x) 

3. 6 = 1 : Z e S5V 

4. 6 = 1 A a = 1/c : Z G TZsv 

There are no nontrivial solutions for general values of a, 6, c. For special 
values we get 

1. a — b/c : Z = Ri iii) 

2. a^b/c : Z ^ Ria{x) 

3. a = c/b : Z E S^v 

4. a = —1, 6 = — c : Z e TJ-sy 

5. a=l, 6 = c : Z e S U {Rl^is)} 
X = ^{a,b) 

There are no nontrivial solutions for general values of a,b. For special 
values we get 

1. a ^ b : Z G Sst 

2. a^-b : Z = i?3.i(-l,-l,l) 

3. a = -6 : Z = {S{x,v) (^t, S{x,y))-^ R^,2{S{x,v) ® S{x,v)) 



where S{x,y) — 




The nontrivial solutions for general values of a are 

1. Z = {A(^ A)Ri,4{x){A(^ A)~^ 

2. Z ^ {A'Si A)Ro,3{A A)-'^ where A 

3. Z = {B(»B)R3,i{x,-x,l){B(»B)-^ 
There are no other solutions for special values of a. 

The nontrivial solutions for general values of 6, c, d are 

1. Z = R2.iic,b/d) 

For special values of b,c,d we get the following solutions {D{q) is 
given by (H)). 

2. b = c = I = -d : Z = D{q)Ri,2{x)D{q)-^ 



72 * 



, where B = 



1 1 
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3. b = d = -c = 1 : Z = D{q)R^%{x)D{q)-^ 
4:. c= -d/b : Z = R2.2(c,x) 

5. c -1 A 6 = -rf = e : Z = D{q)Ro.2D{q)-^ 

6. 6==lAc = rf: Z = {S{x) ® S'(x))i?3.i (d, l/d, l)iS{x) ® S{x))-^ 
wheTeS{x)=(^ ^-J l^d^l 

7. b = c = d=l: Z eSsT 

8. c^l/\b = d^e: Z ^ i?o.i 

• X = A0(6,c,rf,.g) 

There are no nontrivial solutions for general values of parameters. For 
special values we get the following solutions {D{q) is given by (|67|)). 

1. b ~ ed, c = e : Z = i?3.i(e, e, 1) 

2. 6 = ed, c=-e : Z = i?2.i(-e, e) 

3. 6 = -d, c = 1 : Z = D{^)Ri.2{-l + x{\ - 9))D{^)-^ 

4. 6 = d, c = -1 : Z = D(^)Rf,^{~l - x(l - b^))D{^)-^ 

5. b = c = d = I, : Z (z Sst 

6. b ^ 1, c = d = -I : Z = Z ^ {A{x, y) (g) A{x, y))-^ Ro.2{A{x, y) (g) 

where A{x, y) = ^ ^ ^ 

7. d= -b= -c=l : Z = R0.2 

8. c= -b= -d=l : Z = i?o.i 

• X ^ A|^(a,6,c,d), c ^ 

There are no nontrivial solutions for general values of parameters. For 
special values we get 

1. a = c : Z = {A(S) A)R3,i{x,x,l){A(^ A)-'^ 

2. a = -c: Z ^ {A(E) A)Ri,4(x){A(g) A)~^ 

where A = f ^ ^ 

\ —I I 

3. a = c, d^b : Z ^{B{y)®B{y))Rz,i{x,x,l){B{y)®B{y))~^ 
where B(,)=(^(l + ^^) 2/^6(1 - .V2) 

4. a c, d = : Z ={C® C)Ri i{x){C ® C)-^ 

whereC=('^^^(^+^^) 

-eit2 l + e2v2 
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where D(x) 



where E 



5. a = c, d = b : Z = {D{x) ® D{x))RiA{l)iD{x) ® D{x))-^ 
xVb \fh 

eV2 - X exV2 - 1 

6. a = c,d^b: Z = {E <?) E)Ria{x){E ® E)''^ 
iVbei i\fibt2 

1 -eie2\/^ 

7. a = c, d = 6 : Z= [F{y) (g> F(y))i?i.4(l)(F(y) (g> F(y))^i 
V5(eiV2 + ye2) \/& 



where ^(y) = ^ ^ ^ ^ 

a = -c, d ^ -b : Z ^ (G(y) ® G(y))i?3.i(a;,-x,l)(G(y) 

where G{y) ^ ( "^/^^ 



1 y 

9. a = -c, d= -b : Z ^ {H ® H)Ro,3{H (g) H)-'^ 
i\/b ieVb 
1 -e 



where 77 = 



• X^J<J^ia,b) 

The nontrivial solutions for general values of a, 5 are 

1. Z ^ {A(E) A)R2.2{b/a,b/a){A(g) 

2. Z = {A(g) A)Ri.i{a/b){A(g) A)-^ where = 



1 1 

1 -1 



3. Z ={B®B)Ri.2{b/a){B®B)-^ 
where B = \ \ ^ |,(7^ = 1 — bla 

4. Z = (G® C)i?i.2(-&/a)(C'»C')~^ 
where G = ^ J^^ ^ ^ , = 1 + 6/a 

5. Z = {D® D)Rl2{~alb){D® DY^ 
where £> = J ^ , = _i _ a/fo 

6. Z ={E®E)Rl2{a/b){E®E)-^ 

q 
q 



where £; = ( j , g2 ^ -l + a/6 



7. Z = {F® F)Ro,3{F (g) where F = 

There are no other solutions for special values of a, b. 

There is no nontrivial solution of the equations (p 
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The nontrivial solutions for general values of a, r ^ 0, ±1 are 

\. Z = (A(»A)i?i.i(ea/c)(A®^) whereA = (J j , = -^-^ 



• X XJ^ia) 

The nontrivial solutions for general values of a are 

1. Z = {A^ A)R3,i{x,x,l){A(E) A)-'^ 

where ^ = ^ i ^ ^ ' 9^ = ^ 1 

For special values of a we get the following solutions 

2. a = ±l : Z = Ra,i 

3. a = : Z = {B ® B)R-i,i{x,x,l)(B ® B)-^ 

4. a = : Z = (B(g) B)i?i.4(a;)(B(g) where B = 

5. a = : Z = (C(^C)i?i.4(a;)(C0C)-i where C = 



-1 1 
1 1 

—i i 
1 1 



• X = ^(a,6,c,d) 

There are no nontrivial solutions for general values of parameters. For 
special values we get 

1. a — d = : Ze S^v 

2. c ^ b{l + a - d) : Z = Ri.sid - a - 1) 

. X^XJ^{a,b) 

There is no nontrivial solution of the equations (|6 

. X - ^0(6) 

The nontrivial solutions for 6 7^ 0, e are 

1. Z = {A(»A)Ri,4{x){A(»A)~^, 



where A = [ 1 g2 ^ 1 - fe2 4 ^ (6 - iqf 

\ q pq J ' ^ 

For special values of b we get 

2. 6 = 6 : Z = (g) B{x))RQ,2iB{x) (g> B{x))-^ 
where B{x) = 4^ q 

3. 6 = see X^{a = 1) 
X = Ag(a,5, c,e) 

There are no nontrivial solutions for general values of parameters. For 
special values we get 
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1. e = 1,6 = : Z eS^v 

2. e = -1, a = ±1 : Z = Ri,i{i) 

3. e = —1, a = ±1 : Z = Ri,4{x) 

4. e = —1, c — ab : Z G S^v 

5. e = —1, a = —1, c ~ —b : Z G TZsv 

6. e = -1, a=l, c = b: Z eSU 

• X = ^{a,b,c) 

There are no nontrivial solutions for general values of parameters. For 
special values we get 

1. 6 = 0: Z = i?3.i(-l,-l,l) 

2. 6 = : Z= {A{x, y) ® A{x, y))-^ Ro.2{A{x, y) ® A{x, y)) 
where A{x, y) = ( ^ ? 

3. a = c6 : Z G Sst 

There are no nontrivial solutions for general values of parameters. For 
special values we get 



1.6 = 0: Z = {S S)Ro.3{S S)-'^ where S 



1 Vi 

1 -Vi 



2. a = 6 = 0: Z = {S S)R3.i{x, -x, 1){S ® S)-'^ where S 
1 1 



-1 1 ^ 

3. a = 6 = 0: Z ^ (8(8)3) Ri i(x)(S(dS)-^ where S = ( ^ ^ 

V -1 v« 

X = ^(a, 6, c, d, e, f, g, h) and X = Xg(a, 6, c, d, e, /, h). 

Solving the system (|6^) for these two matrices is rather difficult. We can, 

however, simplify the task by the following way: 

Both the matrices solve the equation (^) only for = 1, so that any Z 
that solve (||) for gives solution Z) of the WXZ system. 

This solution is equivalent to [Z^ , X+, 1) due to (p^). It is always possible 
to transform Z^ by WXZ symmetries to the form contained in S. It 
means that for Z ^ 1 we can find the solution of the WXZ by symmetries 
investigating W ^1. It remains to investigate the case Z = 1 or, in other 
words, to solve the equation [X, X, 1] = where X is of the form AO or 
It turns out that the resulting matrices can be always transformea by 
the symmetry transformation X' = {T®1)X{T® 1)"^ to Xi, X2, X3, X4. 
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3.4 List of solutions of the WXZ system 



Combining the results of the subsections 3.2 and 3.3 we get the list of solutions of 
the WXZ system (^-(ph presented below. By construction the list is complete 
up to the symmetries (pl)-(|6l) and (|l|) , (H) , (||) . 

Let us remind that a, b, c, d,g,h ^ C, e, ei, 62, 63 — ±1, and 5, iSgy , iSsy, Sst 
are defined by (|2|)-@. 



Generic solutions: 

1. x = i, z s yj {Rl^is)} 

2. W&S8V, X = 1,1,-1), Z C Ssv {Rl.2{s)} 



3. W ^Sav, X 



4. WeSsv, X = 



/I 













1 

















1 


V 





1 


/ 


/ 1 








\ 





1 

















1 


V 





1 


/ 



Z = i?3.i(a;,a;, 1) 



Z = Ri,i(x) 



^.W€S^v, X ^ diag{l, 1,1, a), Z £ S5V or Z = P 
6. W eS5v, X = 



7. W e Sst, X = 



/ 1 











\ 





1 
















1 







V 





a 


1 




/ 1 











\ 





1 
















1 







V « 








1 





Z e Sst 01 Z = P 



Z e Sst or Z = P 



Non— generic solutions: 



W^ = i?3.i(s,s,l), X 



( 1 








M 











1 








h 





^0 


1 





; 



Z = P 



19 



9. W = R3.i{s,s,l), X 



10. W = R3.i{s,s,l), X = 



11. W = R3.i{s,s,l), X 



12. W = R3.i{s,s,l), X 



13. W = R3.i{s,s,l), X = 



14. W = R3.i{s,s,l), X = 



15. W = R3.iis,s,l), X 



16. W = R3.i{s,s,l), X 



17. W = R3.i{s,s,l), X = 



I 


1 




1 


1 




1 




1 


-1 




\ 




_x 


X 


\ 


-1 




1 


1 


/ 


1 




1 


1 




1 




1 


-1 




1 

L 




_1 


1 


V 


-1 




1 


1 


/ 


1 





















1 










-1 





V 





1 





0/ 


/ 


1 





















1 










-1 





v 





1 








J 


/ 








1 


\ 















a 






1 













V 





1 





) 





Z = R3.i{x,x, 1) 



Z = Ria{x) 



f 10^ 

a 

10 

\ 1/a / 

/ 1 \ 

a 

10 

V 1/a / 

/ 1 \ 

11 

10 

\ c 1 / 

/ 1 \ 
11 

10 

y -1 1 y 



Z = Ri.i{i) 



Z = RiAix) 



Z = P 



Z = Ri.i{i) 



Z = Ria{x) 



Z = P 



, Z = i?3.l(-l,-l,l) 
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18. W = R3.i{s,s,l), X 



19. W = R3.i{s,-s,l), X 



20. W = R3.i{s,-s,l), X 



21. W = R3.i{s,-s,l), X : 



22. W = R3.iis,-s,l), X = 



23. VF = i?3.i(ei,ei,l), X = 









X 












1 



















-1 


X 





/ 




/ 





1 


^ 













-1 




n 


Qj 


n 


n 




V 1 








) 












1 ^ 










1 







n 


Q 


n 


n 




\ -1 








1 




( 








1 ^ 


\ 








ie 







n 


fcCy 


n 







\ 1 




















1 ^ 


\ 








ie 










itq 








1 


V 1 








0, 




( 1 











\ 





ei 


1 











£2 





V 9 








eie2 / 


1 








Z = Ro.2 



Z = P 



Z = R3.i{x, -X, 1) 



Z = Ria{x) 



Z = Ro.3 



Z = -R3.i(e2,e2, 1) 



24. W = R3.iie,e,l), X 



25. W = R3.i{e,€,l), X 



26. W = i?3.i(-l,-l,l), X 



e 1 
10 
\ 9 e J 



Z = R 



•0.1 



/ 1 Q \ 

e 1 

0-10 

V 5 -e y 

f ab c \ 

o -c 

d 6 

V -d 1 / 



, Z — Ro.2 



, Z = P 
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27. W^ = i?3.i(l,l,l), X 



28. W^ = i?3.i (1,1,1), X = 



29. W = R2.i{r,s), X 



30. W = R2.i{r,s), X 



31. W = R2.i{-l/s,s), X = 



32. W = R2.i{e,-e), X 



33. = i?2.i(s,s), X 



34. W^ = i?2.i(«ei,iei), X 





1 













a 


1 


b 












d 





v 


9 


c 


h 


d 


/ 


1 
















a 





6 









d 





V 





.9 





h 


1 


































a 





1 










so 


1 


































a 





1 










sa 




/ 


1 















—s 















a 




V 


1 








/ 1 
















e 
















1 
i 


U 


V 1 








— e 


1 




















1 













a 





9 










sa 




' 1 










- 


-ie\ 


1 










£2 




^ 9 









Z = P 



Z = P 



Z = P 



Z = i?2.i(a,(rsa)-i) 



Z = R2.2{a,x) 



Z = Ri.2{x) 



Z = P 



35. W = R2.iiie,ie), X = 



(I ^ 

Q -ie q Q 

10 

\ a y 



■Z^ = -R2.i(e2, —£2) 



2 = i?1.2(-l+^) 
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36. VF = i?2.i(l,-l), X 



37. W^ = i?2.i(l,-1), X 



38. W^ = i?2.i(l,-1), X 



a 














a 








1 








b 





-1 


b 





a 














a 











1 


6 





1 








-6 


a 














a 











1 


b 





1 








-b 



Z = P 



Z = R2.2{b/a,b/a) 



Z = Ri,i{a/b) 



39. W = R2.i{e,-e), X = 








\ 

ea 

b 



V J-^ -e6 



Z = Ri.2{b/a) 



40. W^ = i?2.i(l,-1), X = 



41. VF = i?2.2(r,s), X = 



42. W^ = ii2.2(r,s), X = 



43. W = R2.2{e,s), X 



44. VF = i?2.2(s,s), X 





( 


a 















a 












h 






\ 


Vi 








/ 1 




































a 







V 


1 





- 


-ar 


/ 1 




































a 







V 


1 





- 


-ar 


/ 1 








] 









e 














1 







1 


V 1 





— e 


) 




( 1 


















1 













a 





\9 








-a/ 



Z = Ro.3 



Z = P 



Z = R2.2{a,x) 



Z = Ri.2{x) 



Z = P 
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45. W = R2.2{s,s), X 

46. W = Ri.i{r), X = 

47. W = Ri.i{r), X = 

48. W = Ri.i{r), X = 
Z = P 

49. W = Ri.i{r), X = 



f 1 i) \ 

s-^ q 

10 

\g -l/s J 

f I \ 

r 1 

-e 

\ e er / 



^ = i?i.2(^-l) 



, Z = P 



(I 

r g 

10 

\ -q -r J 



, Z = Ri.2{^ 1) 



ec 

e 
V -e c 



ea^ 

r— 1 



a 

r+l 
r-1 



\ 

-1 
, r+l 
'r-1 

a / 



c2 = a2 _ 



50. 1/F = i?i.i(i), X 



51. W = Ri.i{i), X = 



52. VF = iii.i(i), X = 



53. W = Ri.i{i), X = 





r+l 
r-1 




















a 


r+l 
r-1 













e 


V 


e 











i?l.l( 


sa/c) 










/ a 


1 - 









1 


—a 
















1 





V 












1 / 


/ e 











\ 




1 


— e 


















1 









V 








-1 


) 




/ 1 








^ 









-1 





















1 




z 


V 





1 


J 






/o 


—i 





M 






i 













z 











1 






V 





1 


/ 







1 \ 

1 

(«-c^)y 



Z = P 



Z = R 



0.1 



Z = Ria{x) 



Z = Ria{x) 
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54. W = Ri.2{s), X = 

55. W = Ri.2{s), X = 

56. W = Ri.2{s), X = 

57. W = i?i.2(s), X = 

58. W = Ri.2{s), X = 

59. W = Ri.2{s), X = 

60. W = Ri.2{s), X = 

61. W = Ri.2{s), X = 
Z = Ro.3 

62. W = i?i.3(u), X = 



/ 1 \ 

10 

e 

V 1 -e y 

/ 1 \ 

10 

10 

V g -1 y 

/ 1 \ 

10 

g 1 



Z = P 



Z = Ri.2{x) 



Z = Rl2{x) 



\ 








-1 / 




( 


1 




















s- 1 

















e 









V 


e 








—es 


) 




/ 


1 













\ 




n 




a{s — 


1) 













1 







V 


a 












/ 


/ 










a(l - 


s) > 









1 

















a 











V 











-1 






/ 


a 










M 









a 














1 










1 




v 





-1 


1? - s 


- 1 


0/ 




/ 


a 

















a 














1/Vi \/ — s 


- 1 


\ 


V* 














/ 


1 










^ 






a 


1 














b 







1 





1 


V 





h — u 


- 1 


—ua 


1 ) 
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Z = P 



z = Ri.2{i 



, z = p 








-Va2-s-i y 



Z = P 



63. W = Ri.3{u), X 



64. W = Ri.3{l), X = 



65. W = RiA{t), X 



66. W = RiA{t), X 



67. W = RiA{t), X = 



68. W = RiA{t), X = 



69. W = RiA{t), X 



70. W = RiA{t), X = 



71. W = RiA{t), X = 



/I \ 

-1/6 1 

b 10 

\ b-u-1 u/b 1 J 

/I ^ 

a 110 

a+b+1 1 

\ a + b-l b 1 J 



Z = R,.3{^-1) 



Z = P 



1 ^ 

a 

e 

\ a y 

f 1 \ 

a 
10 

\ -a / 

^ 10^ 

i 
10 

\ -i / 

/ a 1 - a^ \ 

1 -a 
10 

\ -1 y 



Z = P 



Z = P 



Z = Ria{x) 



Z = P 



f 


e 











\ 




1 


— e 


















1 







V 











-1 


/ 


/ 


1 








\ 









-1 





















1 




V 








1 


) 




/ 


1 








\ 







-1 





















i 




\ 








—i 


0/ 



Z^R 



0.1 



Z = Ria{x) 



Z = Ria{x) 
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72. W = Ri.4{t), X 



73. W = RiAit), X = 



74. W = RiA{t), X . 



75. W = Ria{x), X 



/ 







1 













CI 




1 





\ 1- 


a" 


—a 





/ 








p 








1/p 








e/p 








V 











/ 





1 











a; 


-1 


1 











V -X 


-1 













1 


\ 



ie 
ie 
V 1 / 



76. W = Ria{x), X = 



77. W = Ro.i, X = 



78. W = Ro.i, X 



79. W = Ro.i, X = 



80. W = Ro.2, X = 



/ 1 \ 

ie 

ie 

\ 1 / 

/ 1 \ 

a 

6 e 

\ c -a J 



\ a c 



c 
\ a c 

1 










\ 










-1 







-h 


-1 


/ 








\ 










-1 










-1 


/ 








\ 










e 










— e 


/ 



Z = Ria{x) 



Ro:. 



Z = Ria{x) 



Z = Ro.3 



Z = P 



Z = P 



Ro.: 



Z = P 
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81. W = Rq.2, 



82. W = Ro.2, X 



83. iy = i?o.3, X = 



84. W = Ro.3, X = 



( 


1 




















-1 














n 


n 


_X 


n 

VJ 






V 


q 








1 


} 




/ 





eh 





M 






e 















1 


n 
u 


n 
u 


6 




V 





-1 


1 






/ 


1 








M 









a 


1 















-1 







V 


i 








a ) 




/ 


a 




1 


\ 




b 


a 







1 







i 


-b 




a 


V 


—i 





—a 




b 


/ 


/ 


a 








^ \ 









a 


1 












1 


—a 







V 


i 








ay 





Z = i?, 



0.2 



Z 



Z = P 



Z = P 



85. W = Ro.3., X^ XT ^ X , Z = R 

86. = i?o.4 = P, X = arbitrary matrix, Z = P 



0.3 



4 Conclusion 

It is possible to solve completely the system of equations in the dimension 

two with the assistence of the computer programs for symbolic manipulations. 
Very important tool for both solving and classification of solutions is the use of 
symmetries. 

The invertible solutions of the WXZ system are classified in this paper and 
the authors believe that that they were careful enough so that the solution is 
complete up to the symmetries of the system (|l^)-(^). The rigorous proof of 
the completeness of the solution set is rather difficult task. It is not known for 
pure Yang-Baxter equation either to the best knowledge of the authors. 

The number of the solutions is rather large so that it is not possible to 
investigate in detail all the quantum doubles generated by the Vladimirov's 
procedure. It is therefore necessary to formulate supplementary conditions for 
the solutions that guarantee some suitable properties of the quantum doubles. 
The minimal condition that can be imposed for this goal is that both W and Z 
are second invertible as well |15l. 
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The results of the subsection 3.2, can be reinterpreted as two-dimensional 
representations of the algebras 



W12T1T2 = T2T1W12 
for any invertible solution W of the Yang-Baxter equation. 



(68) 



5 Appendices 

5.1 List of invertible solutions of the Yang— Baxter equa- 
tion 

The important starting point for solving the WXZ system is the knowledge 
of complete set of invertible solutions of the constant Yang-Baxter equation 
[R, R, R] = 0. Up to the symmetries of the solution set of the Yang-Baxter 
equation 

R' ^ piT(g)T)R{T(g>T)-\, R' = R+, R' = R^ , (69) 

it consists ^, |^ of the following (intersecting) subsets parametrized by com- 
plex numbers r,s,t,u,v, 

Solution with three parameters 



R3.i{r,s,t) 



/ 1 \ 

r 

s 

V t y 



rst ^ 



Solutions with two parameters 

/ 1 



R2.i{r,s) 







i?2.2(r, S) 
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/ 1 
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\ 
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1 / 
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s 
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^2.3 (W.^^) = 



/ 1 \ 

110 

M 1 

\ V u 1 1 J 

Solutions with one parameter 
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Ri.2{s) 



Ri.3{u) = 
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RiA{t) = 



Solutions without parameter 



Rf).i = 



Ro.: 



/ 1 \ 

10 

10 

\ 1 1 / 

/ 1 
0-100 

0-10 

\ 1 1 
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-Ro.3 — 



/ 1 i \ 
110 
1-10 
\ i 1 / 
/ 1 \ 
10 
10 
\ 1 / 



5.2 Equivalence cleisses of the 4x4 matrices with respect 
to the transformation X' = (1 (g) S)X{1 S)'^ 

Lemma: Any 4x4 matrix A can be transformed into one of following 14 
forms using transformation 

A' = A(l O S)A{1 O 5)-^ (70) 

where A e C, 1 = diag{l, 1), S e gC{2, C). 
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/ al 0,2 0,3 \ 

bl h2 1 0,3 

cl c3 c4 

\ dl cl d3 d4 I 

I al a2 a3 \ 

61 62 a3 

a + 1 c3 c4 
\ a - 1 d3 d4 y 
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All 
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Proof: In order to prove previous lemma it is convenient to write matrix A 
divided into square blocks, i.e. 



A = 



All Ai2 
All A22 



(71) 



XSAiiS ^ under tranformation 



Blocks All, • • ■ , A22 are transforming like A[j — ^v^^ij^^ 

([70I). It is possible to convert one chosen block into Jordan's canonical form, 
i.e. one of the following forms: 



1) 



a 
a 



2) 



a 
1 a 



3) 



a + 1 






a-1 



(72) 



( see Jordan's theorem in linear algebra ). Case 1) is of course invariant w.r.t. 
([to]), in case 2) it is possible to convert another arbitrary block into upper 
triangular ( with non-diagonal element equal 1 ) or lower triangular ( with 
identical diagonal elements ) or diagonal block. Using these properties, one 
can find given list of classes, firstly simplifing upper right block, secondly using 
remaining symmetries to simplify lower left block, then lower right block and 
finally upper left block. 
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5.3 List of X matrices 

Solutions of the equation (0) can be written in terms of the following matrices. 



Xi{a,b,c) = 



X2{a,b,c) = 



X3{a,b,c) 



^4(0, 6, c) 
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Xio{b, c, d,g) 
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X2oia,b,c) 



X2i{a,b) 



X22{a,b, c, d,e,f,g,h) 



X23{a,b, c, d,e,f,g,h) 
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